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We evaluate the matrix elements of the four-quark operators which contribute to the lifetimes of B-mesons and 
the Af,-baryon. We find that the spectator effects are not responsible for the discrepancy between the theoretical 
prediction and experimental measurement of the ratio of lifetimes t(Ai,)/t(B). 
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1. INTRODUCTION 

Inclusive decays of heavy hadrons can be stud- 
ied in the framework of the heavy quark expan- 
sion, in which, for example, lifetimes are com- 
puted as series in inverse powers of the mass of 
the &-quark fjj]. For an arbitrary hadron H 



-\H) 



where 



G 2 F ml \V cb \ 2 ^ _, 
2^ c * m b 



1927T 3 2m H 



(1) 



j>0 



• Co corresponds to the decay of a free-quark 
and is universal. 

• ci is zero because the operators of dimen- 
sion four can be eliminated using the equa- 
tions of motion. 

• C2 can be estimated and is found to be 
small. 

• C3 contains a contribution proportional to 



(H\ bTqqTb \H) 



(2) 



and is therefore the first term in the expan- 
sion to which the interaction between the 
heavy and the light quark(s) contribute. Al- 
though this is an 0(m^ 3 ) correction, it may 
be significant since it contains a phase-space 
enhancement. 

The aim of our lattice simulation is to compute 
C3 for B-mesons and the A^-baryon, in order to 



check whether spectator effects contribute signif- 
icantly to the ratios of lifetimes for which the ex- 
perimental values are: 



t(B~ 



t(B°) 
t(B°) 



= 1.06 ±0.04 



= 0.78 ± 0.04 . 



(3) 
(4) 



The discrepancy between the experimental 
value in eq. (f|) and the theoretical prediction 
of r(A 6 )/r( J B°) = 0.98 (based on the Operator 
Product Expansion in eq. (|l|) including terms in 
the sum up to those of 0(m^ )) is a major puzzle. 
It is therefore particularly important to compute 
the 0(m^~ 3 ) spectator contributions to this ratio. 



The ratios in eqs. (||) and 
in terms of 6 matrix elements: 

t(B- 



can be expressed 



t(B°) 

r(B°) 
where 2 

B\ = 
B 2 EE 
e\ EE 



a + axSi + a 2 e 2 



= b + b\E\ + b 2 e 2 



b 3 Li 4 



a 3 B 2 (5) 
(6) 



b A L 2 



8 (B\ b^Lqq^Lb \B) 
f%m B 2m B 

8 (B\ bLq-qRb \B) 
J^Tb 2^ 

8 (B\ ~try> 1 Lt a qqj t ,Lt a b \B) 
f%m B 2^ 



B 



(7) 
(8) 
(9) 



2 In terms of the parameters B and r introduced in ref. [|2 
r = —6Li 



B 



-2£ a /£i - 1/3 
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Li 



Ln = 



8 (B\ bLt a qqRt a b \B) 
f B m B 2m B 

8 (A| bj'Lg q lf2 Lb |A) 
f B m B 2m A 

8 (A| Tr/»Lt a qq>y fi Lt a b |A) 



(10) 



quark propagators are computed using a stochas- 
tic method Ik) based on the relation 



f B m B 2m A 
and the coefficients and are given by 





value 




value 




+ 1.00 


bo 


+0.98 


ai 


-0.697 


h 


-0.173 


0,2 


+0.195 




+0.195 


d3 


+0.020 


b 3 


+0.030 


CL4 


+0.004 




-0.252 



The values in the table correspond to opera- 
tors renormalized in a continuum renormalization 
scheme at the scale [i — tub- Their matrix ele- 
ments are obtained from those in the lattice reg- 
ularization by perturbative matching (see the ap- 
pendix) . 

2. B DECAY 

The matrix elements Bi,B 2 ,£i,e2 are com- 
puted on a 24 3 x 48 lattice at (3 — 6.2 (correspond- 
ing to a lattice spacing a -1 = 2.9(1) GeV) using 
the tree-level improved SW action for three values 
of k = 0.14144, 0.14226, 0.14262 and are then ex- 
trapolated to the chiral limit (k c = 0.14315) |J. 
We find 

B 1 = 1.06 + 0.08 (13) 
B 2 = 1.01 + 0.06 (14) 
Ei = -0.01 + 0.03 (15) 
e 2 = -0.02 + 0.02 (16) 

which implies that 

T<yB } = 1.03 + 0.02 + 0.03 (17) 
t(B v ) 

in agreement with the experimental value (|j|). 

3. A DECAY 

The computation the baryonic matrix elements 
L\ and L 2 is a little more difficult. We have per- 
formed an exploratory study in which the light 



(11) My 1 = / [#](M,- fc <fe)*&e-^ M+M ^ (18) 



(12) 



(rather then using "extended" propagators). 

The matrix elements are computed on a 12 3 x 
24 lattice at /3 = 5.7 (corresponding to a lattice 
spacing a -1 = 1.10(1) GeV) for two values of k. 
We therefore do not attempt an extrapolation to 
the chiral limit (k c = 0.14351) but present results 
seperately for each value of k. We find: 

-0.30 + 0.03 (k= 0.13843) 
-0.22 + 0.03 (k= 0.14077) 



Li 



(19) 



0.23 + 0.02 (k = 0.13843) 
0.17 + 0.02 (k = 0.14077) 



(20) 



(21) 



which implies that (neglecting the systematic er- 
ror due to the chiral extrapolation) 

r(A b ) _ f 0.91 + 0.01 (k= 0.13843) 
t{B°) ~ \ 0.93 ± 0.01 (k = 0.14077) 

These results imply that spectator effects are not 
sufficiently large to explain the discrepancy be- 
tween the theoretical prediction and experimental 
result in eq. (Q). 

4. CONCLUDING REMARKS 

We find that the matrix elements of the 4-quark 
operators (p^-^6[) satisfy the vacuum saturation 
hypothesis remarkably well. A similar feature is 
true for the AB = 2 operators which contribute 
to B—B mixing. We do not have a good under- 
standing yet of this phenomenon. 

The results for the mesonic matrix elements 
lead to a prediction for the ratio of the lifetimes 
of the charged and neutral mesons which is in 
agreement with the experimental result in eq. (g) . 
Our study of the baryonic matrix elements in- 
dicates that spectator effects are not sufficiently 
large to explain the experimental ratio of life- 
times in eq. (^). This discrepancy between the 
theoretical prediction and the experimental mea- 
surement remains an important problem to solve. 
We do stress, however, that our calculations are 
exploratory, and a more precise simulation is nec- 
essary, in particular to allow for a reliable extrap- 
olation to the chiral limit. 
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Table 1 

Lattice perturbative coefficients and operators 



Pi 


Pi 


Qi 


Qi 


-f log(AV) + 8.46 


bTqqTb 


+ 1 


bTqbVq 


-log(A 2 a 2 ) + 6.19 


bt a Tt a qqTb + bTqqt a Tt a b 


+1 


bt a Tt a qbfq + bTqbt a Yt a q 


-6.89 


bt a j°T'y t a qqfb + bTqqt a j°f j a t a b 


+1 


bt a ^ Tj°t a qbfq + bTqbt a j a fj°t a q 


21og(A 2 a 2 ) -4.53 


bt a Tqqft a b 


-1 


bt a Tqbt a fq 


log(A 2 a 2 ) - 6.19 


bTt a qqft a b + bt a Tqqt a Tb 


-1 


bTt a qbt a fq + bt a Tqbft a q 


-6.89 


bT^°t a qqfj t a b + bt a j°Tqqt a -f fb 


+1 


bTj a t a qbt a j°fq + bt a j°Tqbfj°t a q 


-log(A 2 a 2 ) + 5.12 


bTt a qq-y°fb 


-1 


bTt a qbTt a q 


-ilog(A 2 a 2 ) + 1.05 


bTa^t a qqt a cr"^Tb 


-1 


bTa^ u t a qbTc7^ u t a q 


-2.43 


bT^t a qqt a ^fb 


+ 1 


bT^t a qbY^t a q 



The analytic expression for pi can be found in ||j . 



ACKNOWLEDGEMENTS 

It is a pleasure to thank Chris Michael, Hart- 
mut Wittig, Jonathan Flynn, Luigi Del Debbio, 
Giulia De Divitiis, Vicente Gimenez and Carlotta 
Pittori for many helpful discussions. This work 
was supported by PPARC grants GR/L29927, 
GR/L56329 and GR/I55066. 

APPENDIX: 1-loop lattice perturbative 
corrections to 4-quark operators 

The most difficult component of the evaluation 
of the 1-loop perturbative matching between four- 
quark lattice operators and those renormalized in 
a continuum scheme is the perturbative expan- 
sion on lattice. In this appendix wc present the 
corresponding results for generic operators Po and 
Qo (whose matrix elements contribute to lifetimes 
and B-B mixing respectively): 



P = bTqqTb, (AB = 0) 
Qo = bTqbfq, (AS = 2) 



(22) 
(23) 



r®r represents an arbitrary spinor and color ten- 
sor. These operators mix under renormalization 
with other 4-quark operators, listed in table |l| [| 

a s (a _1 ) 



1 loop 



P, 



Ql loop 




Po 



Qo 



47T 
47T 



i 

^Piq l Qi 



(24) 



(25) 



The Feynman rules correspond to the tree-level 
SW-improved action for massless light quarks and 
with a small gluon-mass (A) as the infrared reg- 
ulator. The dependence on A, of course, cancels 
when the corresponding continuum calculation is 
combined with the lattice one. 
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